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Abstract 

Lattice J\f = 1 super- Yang-Mills theory formulated using Ginsparg- Wilson fermions provides a 
rigorous non-perturbative definition of the continuum theory that requires no fine-tuning as the 
lattice spacing is reduced to zero. Domain wall fermions are one explicit scheme for achieving 
this and using them we have performed large scale Monte Carlo simulations of the theory for 
gauge group SU{2). We have measured the gaugino condensate, static potential, Creutz ratios 
and residual mass for several values of the domain wall separation Lg, four-dimensional lattice 
volume, and two values of the bare gauge coupling. With this data we are able to extrapolate the 
gaugino condensate to the chiral limit, to express it in physical units, and to establish important 
benchmarks for future studies of super- Yang-Mills on the lattice. 
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I. INTRODUCTION 



The lattice formulation of supersymmetric theories necessarily breaks most or all^ of the 
continuum supersymmetry, since the translation group, a subgroup of the super-Poincare 
group, is broken to a discrete subgroup. This is intimately related to the failure of the 



Leibniz rule on the lattice jsj. Hence all relevant and marginal operators that are allowed 
by lattice symmetries will be generated radiatively, so that generically the long distance 
effective theory will have supersymmetry badly broken. 

On the other hand, the only relevant or marginal operator allowed in a gauge invariant 
lattice formulation of pure M = 1 super- Yang-Mills ^\ (SYM) with hypercubic symmetry 
is the gaugino mass term, as was emphasized long ago in the analysis of jsl. By using 
Ginsparg- Wilson fermions j^, the associated lattice chiral symmetry Q protects against 
additive renormalizations of the gaugino mass in the continuum limit. Hence the desired 
continuum theory is obtained without fine-tuning of counterterms, merely by setting the 
bare fermion mass to zero. 

A viable lattice formulation provides nonperturbative information regarding SYM, the 
foundation for all nonabelian four-dimensional supersymmetric gauge theories. Given the 
unique features of supersymmetric theories, it makes sense to study their strong dynamics 
by as many means as possible. It is quite exciting that the lattice approach can be brought 
to bear on A/" = 1 SYM in a rigorous and reliable way, since it would be enlightening to 
study the nonperturbative aspects of this theory in detail using numerical techniques. For 
instance, information on the spectrum and renormalization of nonholomorphic operators 
would be welcome. 

It is essential to quantify the size of lattice artifacts, since these lead to deviations from 
the continuum physics. One of the purposes of the present study is to determine the regime 
of lattice parameters where continuum SYM makes its appearance (to a good approxima- 
tion), and to characterize the difficulty of performing Monte Carlo studies in that limit, 
given the algorithms and computing resources that are currently available. The study that 
we present is representative of what can be achieved with a dedicated world-class supercom- 



^ In cases with extended supersymmetry it is sometimes possible to preserve a nilpotent subalgebra (see for 
example the reviews [l|, Q and references therein). For minimal supersymmetry, such as is studied here, 
all of the generators are broken. 
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puting resource — of order^ 1 Teraflop/s x year (30 Million IBM BlueGene/L core hours) on 
Rensselaer's Computational Center for Nanotechnology Innovation — and highly optimized 
parallel code (a modified version of the Columbia Physics System). It will be seen that we 
are able to obtain reliable extrapolations to the chiral limit (vanishing gaugino mass), but 
that the study was quite demanding and could not have been performed with anything less 
than the resources just described. 

Whereas much is known about the vacuum structure of SYM by continuum methods 

0, y, Q, Q, Q, Q, 

nothing is known about its long distance dynamics and spectrum. 
Furthermore, lattice methods have the ability to reveal far more detail about the vacuum, 
as one can subject it to a configuration by configuration analysis, as has been done for pure 
Yang-Mills with considerable success. Finally, we are interested in the effects of a small, 
nonzero gaugino mass in the nonperturbative regime. The explicit, but controllable, chiral 
symmetry breaking of the fermion discretization we use here allows us to explore the impact 
of a small gaugino "soft mass" on nonperturbative quantities such as the condensate, the 
string tension and the shape of the static potential. This is in the spirit of a large number of 

nnfinnnr 

proposals made over the last several years PJ, ll5|, M, UJl ll8|, UJ, |20l . Ultimately, the impact 
of nonzero gaugino mass on the spectrum of bound states will emerge from high statistics 
studies that are beyond the scope of the present work. It will then be possible to compare 
the lattice data to the references just cited. 



One implementation of Ginsparg- Wilson fermions is domain wall fermions (DWF) 21 
22l ] in the limit of infinite separation between the walls, — oo. Besides the absence 
of nonperturbative fine-tuning of the gaugino mass, DWF have the advantage that the 
fermion measure is real, positive and the square root of the determinant which enforces 



the 



23 



Vlajorana condition is analytic with a phase that is independent of the gauge fields 



24| . These three features are all lacking in the Wilson fermion formulation that was 



applied in the only concerted effort to date to study SYM on the lattice by the DESY 



Miinster-Roma collaboration 



25 



28 



301] and to a lesser extent Donini et al. 31 1 



(Recently, this program has been revived [32].) Our research, which has already appeared 
in preliminary form [33] , is in some sense a continuation of the work of Fleming, Kogut and 



^ This is the actual computing power brought to bear. The theoretical computing power utihzed was ten 
times this. 
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Vranas (FKV) 3J] who pioneered the use of DWF for studying TV = 1 SYM. Similar work 



has been initiated by Endres 



351]. What sets the present study apart is that an extensive scan 



of the domain wall separation Lg and measurement of the residual chiral symmetry breaking 
mass m^cs was done at two values of the bare lattice gauge coupling {(3 = A/g'^ = 2.3 and 
2.4) and spatial/temporal volumes (L^ = 8^ and 16^; T = 16,32). This has allowed us to 
obtain a reliable chiral extrapolation (mres 0), and a preliminary view on what occurs as 
we take the continuum, theormo dynamic limit {l3,L,T —>■ oo). 

The lattice formulation that is used in this study has already been described by FKV |3J] ; 
it is reviewed in Appendix El In brief, the lattice employs Shamir DWF 22|] in the adjoint 
representation of SU(2), and the one-plaquette Wilson gauge action. The Majorana con- 
dition is imposed through a square root on the fermion determinant, which as mentioned 
above is analytic and introduces no gauge field dependent sign ambiguity 



23 



24|. 



All results reported in this article utilize a domain wall height mo=1.9, as in the FKV 
simulations, though we will comment briefly on some tests we did at other values of mg. 
Lattice configurations were generated with a dynamical fermion mass mj = 0, so that 
the finite size of the fifth dimension, parameterized by Ls, was the sole infrared regulator, 
through the corresponding additive mass correction m^ps (reviewed in Appendix [B]) , which 
is a measure of residual chiral symmetry breaking 36] . One does not expect an optimal 
value of mo to exist 3^], but for stronger couplings, the range is rather narrow. Our mres 
measurements below show that we are in the correct phase — the explicit chiral symmetry 
breaking decreases as Lg is increased — for mo = 1.9. Finally, using the lattice configurations 
that we generated, we computed mres for DWF propagators with other values of mo. We 
found that mres could be lowered slightly by increasing mo toward the critical value 2, and 
that decreasing mo from 1.9 increased mres- The decrease in mres by increasing mo was not 
significant, so we did not pursue the issue further. 

In Section [III we give our results for the bare gaugino condensate for various couplings 
(3 = 4/(7^ and domain wall separations Lg. Then in Section IITII we summarize our findings 
for "gluonic" observables (i.e., those related to the nonabelian gauge bosons), principally 
the string tension and thereby the Sommer parameter ro/a. Next in Section [IVl we discuss 
the chiral extrapolations of the gaugino condensate based on our simulation results. We 
conclude in Section |Vl In addition to the two appendices mentioned above. Appendix O 
describes technical aspects of the simulation that may be of interest. 
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/3 


2.1 


2.2 


2.3 


2.4 


2.5 


(AA)a3 


0.007494(5) 


0.007719(5) 


0.007051(5) 


0.00499(6) 


0.003043(5) 



TABLE I: A scan of the condensate versus f3 for the 16^ x 32 lattice with domain wah separa- 
tion Lo = 16. 



A 



0.000 ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' 

2.0 2.1 2.2 2.3 ^ 2.4 2.5 



FIG. 1: Condensate vs. /3 for 16^ x 32 lattice with Lg = 16. 



II. BARE GAUGING CONDENSATE 



We have validated our simulation code by comparing to FKV at several points. We 
obtained results that agree with FKV, to within 1% statistical errors. 

A summary of all results obtained here for the gaugino condensate (AA) is given in Ta- 
bles [H [TTl and [TTll The residual chiral symmetry breaking is parameterized through mres 36| . 
which we briefly review in Appendix |Bl Measurements were conducted on large and small 
lattice volumes; it can be seen that in lattice units the finite-size dependence is mild or in- 
significant for (3 = 2.3 but quite noticeable for f3 = 2.4. This is sensible, given that f3 = 2.4 
corresponds to a finer lattice spacing, and hence the physical volumes are smaller. Simula- 
tions on 16'^ X 32 volumes with = 48 are in progress and will be presented elsewhere. 

For the Lg = 16 lattices, which are relatively inexpensive, a scan over (3 was performed, 
with results given in Table [T] and shown in Fig. [1] The vanishing extrapolated value at 
/? ~ 2.7 is apparently due to finite-size effects that cause the system to deconfine. 

We have measured the condensate at other values of Lg using a sea-Ls/valence-L^ ap- 
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V xT 


Ls 




(AA)a3 




(AA)r3 


8^ X 8 


16 


0.158(5) 


0.00711(7) 






8^ X 32 


16 


0.181(3) 


0.00703(4) 


0.75(13) 


0.51(27) 


16^ X 32 


16 


0.184(2) 


0.007051(5) 


0.668(10) 


0.337(11) 


8^ X 32 


24 


0.1541(15) 


0.005112(8) 


0.610(97) 


0.32(15) 


16^ X 32 


24 


0.1564(17) 


0.005321(9) 


0.546(55) 


0.226(68) 


8^ X 32 


32 


0.1319(12) 


0.004321(11) 


0.501(69) 


0.24(10) 


16^ X 32 


32 


0.143(2) 


0.00445(2) 


0.483(58) 


0.172(61) 


8^ X 32 


40(1) 


0.1183(54) 


0.00383(3) 






8^ X 16 


48 


0.1043(17) 


0.003563(20) 


0.361(31) 


0.148(37) 


8^ X 32 


48 


0.1071(10) 


0.003551(11) 


0.409(31) 


0.198(45) 


8^ X 32 


64 


0.08864(84) 


0.003164(10) 


0.300(35) 


0.122(42) 



TABLE II: The gaugino condensate (AA) and residual mass rrires for various lattice sizes and Ls 
values, all at /? = 2.3. The Lg = 40 value, with an "(I)" after it, is obtained by the interpolation 
method described in the text, and summarized in Table IIV[ Values in units of the Sommer pa- 
rameter ro are also shown, for those cases where the potential was measured (in particular, for all 
points that are included in the chiral extrapolation fit). The Lg = 16 data was not included in 
the linear chiral extrapolation fit, because these points had too much curvature (with respect to 
mres«) associated with them. 



proach. The condensate was measured using DWF with L™'' on top of dynamical lattices 
produced using a nearby U^^. Performing this for U'^^ values on either side of L™'' yields 
robust interpolated (I) results, as can be seen in Table [IVl These are then used in our fits of 
(AA) vs. mres, together with the strictly dynamical {Uf^' = L^'^^) measurements of Tables HTl 
and lllll Interpolations ("I") are taken halfway between the results. Half the difference plus 
the statistical errors added in quadrature is used as an error estimate for the interpolation. 

We also use the results of Section IIIIBI below to express m-res and (AA) in terms of 
the Sommer scale ro 



38| . Note that the P = 2.4 value of mres'^o at Lg = 48 indicates 
that the effective gaugino mass (which should be approximately equal to mj-es) is roughly 
1/4 the inverse Sommer scale, so that we are beginning to enter the chiral regime where 
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V xT 


Ls 




(AA)a3 




(AA)r3 


8^ X 32 


16 


0.080(2) 


0.004839(15) 


0.547(30) 


1.55(22) 


16^ X 32 


16 


0.0969(8) 


0.00499(6) 


0.5355(66) 


0.842(25) 


8^ X 32 


24 


0.0601(15) 


0.003293(17) 


0.417(26) 


1.10(18) 


16^ X 32 


24 


0.0838(17) 


0.00389(8) 


0.385(35) 


0.38(10) 


16^ X 32 


28(1) 


0.0721(33) 


0.003452(45) 






8^ X 32 


32 


0.0486(12) 


0.00269(2) 


0.296(15) 


0.61(08) 


16^ X 32 


32 


0.0653(15) 


0.003330(12) 


0.313(33) 


0.37(11) 


8^ X 32 


40(1) 


0.0390(24) 


0.00234(8) 






8^ X 32 


48 


0.0328(9) 


0.002165(18) 


0.224(17) 


0.69(15) 



TABLE III: Results similar to Table HIl except that these are for (3 = 2.4. 



/? 


V 


T 








(AA)a3 


2.3 


83 


32 


40 


32 


0.117(4) 


0.003818(9) 


2.3 


83 


32 


40 


48 


0.1196(10) 


0.003843(9) 


2.3 


83 


32 


40 


I 


0.1183(54) 


0.00383(3) 


2.4 


163 


32 


28 


24 


0.0707(13) 


0.003407(3) 


2.4 


163 


32 


28 


32 


0.0734(15) 


0.003496(3) 


2.4 


163 


32 


28 


I 


0.0721(33) 


0.003452(45) 


2.4 


83 


32 


40 


32 


0.0381(10) 


0.002284(13) 


2.4 


83 


32 


40 


48 


0.0398(11) 


0.002397(17) 


2.4 


83 


32 


40 


I 


0.0390(24) 


0.00234(8) 



TABLE IV: The valence/sea results and interpolations ("I"). Here -L™ is the value used for the 
measurements and L^^^ is the value used in the dynamical fermion simulations. 

supersymmetry is well approximated. On the other hand, it can be seen that rrij-esro is 
unpleasantly large for (3 = 2.3 with < 32, and likewise the condensate in physical 
units is small compared to the (3 = 2.4 results. Clearly (3 = 2.3 is further away from the 
supersymmetric limit due to the coarser lattice. On the other hand it can be seen that the 
(3 = 2A data shows a marked volume dependence due to the smaller physical "box" that 
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the states must squeeze into. 

We note that in the present context the mres measurement coming from the midpoint 
"pion" propagator calculation is a quenched probe of explicit chiral symmetry breaking, since 
in the dynamical theory (i.e., the one pertaining to the lattice action that is used to generate 
configurations) we do not have two "valence quarks" and a nonanomalous continuous chiral 
symmetry that would give rise to pseudo-Nambu-Goldstone bosons. Rather, the (adjoint) 
pions that are measured in the midpoint calculation of pseudo-Nambu-Goldstone 
bosons of an SU{A\3) graded chiral symmetry that is spontaneously broken in the chiral 
limit, and explicitly broken at finite Lg, as we now explain. 

The SU{4:) subalgebra of the graded Lie algebra S'f/(4|3) is associated with the gaugino 
(a Majorana fermion) plus three quenched Majorana fermion probes. The total of four 
Majorana fermions with degenerate mass is equivalent to two Dirac fermions in the adjoint 
representation, with a resulting SU{4) chiral flavor symmetry in the chiral limit. The fact 
that three of the four Majorana degrees of freedom are quenched is equivalent to introducing 
three Majorana ghosts, also with the same mass. As in partially quenched QCD, the ghosts 
cancel the contribution of the nondynamical fermions to the functional integration measure. 
Also analogous to quenched QCD, there is a graded Lie algebra that relates the fermion 
and ghost fields, namely S'[/(4|3) in the present case. The PCAC mass (mres) associated 
with this Nf = 2 adjoint-Dirac fermion chiral symmetry breaking is a good probe of the 
DWF chiral limit, for the same reasons that it is a solid tool in quenched DWF-QCD studies. 
Investigations of the effective theory description of the SU{4:\3) algebra, from the theoretical 
perspective as it relates to DWF, are in progress 39|. Finally we note that in recent spectrum 
studies of one-flavor QCD, a similar non-singlet flavor current was utilized with success 40 |. 



III. GLUONIC OBSERVABLES 



One of the interesting features of SYM is that it is a theory with dynamical fermions 
that do not screen static sources in the fundamental representation. In contrast to QCD, 
it has true confinement in the sense of an area law and no string breaking. (Recall that 
the gauge action is expressed in terms of fundamental links, so we are not studying the 
SU{2)/Z2 = 5*0(3) gauge theory, wherein fundamental sources would have an ambiguous 
meaning.) On the lattice, we can therefore study a very interesting static potential — one 
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/9 


x(i,i) 


X(2,2) 


X(3,3) 


X(4,4) 


X(5,5) 


X(6,6) 


2.1 


0.6423(6) 


0.5242(16) 


0.459(7) 


0.4(2) 






2.2 


0.580(2) 


0.427(5) 


0.330(10) 


0.12(5) 






2.3 


0.51051(11) 


0.3091(4) 


0.1997(12) 


0.163(11) 


0.17(7) 




2.4 


0.45966(7) 


0.2346(3) 


0.1227(9) 


0.078(3) 


0.052(14) 


0.19(7) 


2.5 


0.42493(8) 


0.1969(3) 


0.0896(5) 


0.0507(17) 


0.036(7) 





TABLE V: Creutz ratios for the 16^ x 32 x 16 lattice. 



with chiral fermions and a nonvanishing string tension. These features of SYM will be 
presented here, illustrating the special ability of the lattice approach: to conduct detailed 
studies of the nonperturbative aspects of the theory that the continuum methods cannot 
touch upon. 



A. Creutz ratios 



Here we look at Creutz ratios 



4l| as a probe of the string tension in lattice units, a a?, as 
well as to delineate the scaling regime where the continuum limit may be extracted. Results 
for the 16'^ x 32 x 16 lattice are shown in Table |V] and Fig. [21 Although the errors are 
somewhat large, scaling is clearly setting in at around (3 ~ 2.4 as can be seen by the x(4, 4) 
ratios, which lie quite close to the 2-loop curve. Much beyond that (3, finite size effects 
will take over and it is necessary to move to a larger lattice. For this reason, most of our 
simulations have been performed aX (3 = 2.3 and (3 = 2.4. 



B. Static potential 

The static potential was obtained by measuring Wilson loops with one side of length t in 
the temporal direction, according to standard methods. Having obtained V{r)a from fitting 
the exponential decay in time, we next analyze the potential in terms of the distance ra. 
We fit the data to the standard form 

V{r)a = Voa + aa'^{r/a) - (3.1) 
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1.000 
0.500 



0.100 
0.050 



0.010 — 
0.005 — 



2.0 



it 

[] 



2. 1 



2.2 



+ = x(i,i) 

X = X(2,2) 
o = X(3,3) 



2.3 
/3 



i 



2.4 



□ ='x(4,4) 
^ = X(5,5) 
M = x(6,6) 



2.5 



2.6 



FIG. 2: Creutz ratios for the 16^ x 32 x 16 lattice. The dashed hne indicates the 2-loop prediction 
for the dependence a^(/3), apart from an overall normalization that is determined by requiring that 
the curve pass through the data. 



as was also done by the DESY-Miinster-Roma collaboration in 25|, and in the recent 



work 



351]. We obtain the Sommer parameter r^/a from this fit, using the formula 



a 



1.65 — a 



(3.2) 



This approach to the determination of ro/a has some sensitivity to the range of radii that 
is fit, and obviously depends on what form we assume for V{r)a. The force method of 
Sommer |38| improves on this by avoiding a static potential fit, relying instead on tree- 
level improved finite differences. Preliminary studies that we have conducted indicate that 
smearing of gauge links and higher statistics will be indisposable to a successful application 
of the method, and so we will perform this in a future study. 

The results of our static potential fit and the derived quantities are presented in Tables |VT] 
and IVIIi Note that the values of the string tension in physical units, otq, are quite uniform 
for fixed f3, and have small errors. For the L = 16 results, the fits were also done using the 
same set of Wilson loops as in the L = 8 case, denoted "L = 8 method", so that dependence 
choice of Wilson loops could be controlled for, and therefore ruled out as a spurious source 
of finite size dependence. It also gives a clearer picture of the systematic uncertainties 
associated with the measurement and fitting methods for the static potential. 
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These results are to be compared with Table 1 of [35|. There, a nonzero fermion mass 
m/ = 0.02 was used, Ls = 16, and somewhat different fit ranges for r,t were employed. In 
particular, our fits include the points r/a = 1, -\/2, which have very small errors and thus 
strongly influence the fit. Thus for = 16 we have also performed a fit with rmin = as 
was done in [35|], as can be seen in the second = 16 entry for /3 = 2.3, 2.4 in Tables IVTl and 
IVIII Our results for the fit quantities at /3 = 2.3 are in good agreement once this restriction 
is imposed. For our other Lg values we have far fewer samples, as larger Lg comes at greater 
computing cost. The degradation of statistical errors that results if we exclude the r = 1, 
joints is unacceptable, which is why we do not quote results with the as 



35( 1 for the other Lg values. On the other hand, it can be seen from the Lg = 16 results 



that the choice of r^i^ only has a 10% effect on the rg/a estimate, so that the choice of r^in 
is not crucial to the broad picture that we are after in this preliminary work. The (3 = 2.4 
results are also in reasonable agreement with [3^, comparing to the numbers we obtain at 
16^ X 32, Lg = 16, with r^i^ = ^/3, and keeping in mind the nonzero mj in 35|]. 

Above, we have used the results of Tables I VI I and IVIII to scale the residual mass and 
condensate to tq units. (Note that the ro/a values with identical lattice parameters were 
used in this procedure, rather than a uniform rg/a value across all mres and (AA).) With the 



string tension in hand, we now see that the energy scale of confinement a/ot^ 1.4 lies above 
the explicit chiral symmetry breaking scale mres'^o by a factor of 1.8 to 3.8 for (3 = 2.3, and 3.0 
to 5.2 for j3 = 2.4. This is consistent with the observation that the string tension results in 
Tables NH and I VI J are insensitive to the range of Lg values displayed there, when expressed in 
physical units (otq). That is, confinement dynamics are to a good approximation decoupled 
from the explicit chiral symmetry breaking. Since the lowest lying excitations of SYM are 
glueballs and superpartners, the gap associated with confinement should also decouple these 
states from the explicit chiral symmetry breaking. Thus it appears that we are well into 
the regime where the spectrum reflects super symmetry, and it will be quite interesting to 
examine the spectrum in order to check whether or not this is true — something we will 
do in future work. In addition, it gives us confidence that we are performing the chiral 
extrapolation of the condensate (next section) correctly, where the data points are dominated 
by the physics of the supersymmetric theory. 
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V xT 


Ls 


Foci 




a 


ro/a 


arl 

U 


method 


8^ X 32 


16 


0.717(83) 


0.074(28) 


0.368(55) 


4.16(73) 


1.282(55) 


L = 8 


16^ X 32 


16 


0.6533(80) 


0.1004(25) 


0.3271(57) 


3.630(39) 


1.3229(57) 


L = 16 


16^ X 32 


16 


0.489(45) 


0.1367(88) 


0.159(55) 


3.303(50) 


1.491(55) 


L = 16,r > 


8^ X 32 


24 


0.718(89) 


0.082(30) 


0.371(60) 


3.96(63) 


1.279(60) 


L = 8 


16^ X 32 


24 


0.752(70) 


0.102(24) 


0.411(46) 


3.49(35) 


1.239(46) 


L = 16 


16^ X 32 


24 


0.696(67) 


0.119(22) 


0.372(46) 


3.27(24) 


1.278(46) 


L = 8 


8^ X 32 


32 


0.748(82) 


0.087(27) 


0.400(55) 


3.80(52) 


1.250(55) 


L = 8 


16^ X 32 


32 


0.745(90) 


0.109(31) 


0.412(59) 


3.38(40) 


1.238(59) 


L= 16 


16^ X 32 


32 


0.635(50) 


0.146(17) 


0.338(33) 


3.00(14) 


1.312(33) 


L = 8 


8^ X 16 


48 


0.706(68) 


0.107(22) 


0.372(47) 


3.46(29) 


1.278(47) 


L = 8 


8^ X 32 


48 


0.768(47) 


0.085(15) 


0.414(33) 


3.82(29) 


1.236(33) 


L = 8 


8^ X 32 


64 


0.680(94) 


0.113(32) 


0.353(63) 


3.38(39) 


1.297(63) 


L = 8 



TABLE VI: Gluonic observables obtained from the static potential for /? = 2.3. We note that 

there is no sign of volume dependence in the string tension results expressed in physical units oTq. 
Errors estimates are obtained from a jackknife analysis of fits. 



V xT 


Ls 






a 


ro/a 


arl 


method 


8^ X 32 


16 


0.617(11) 


0.0292(30) 


0.2857(83) 


6.84(33) 


1.3643(83) 


L = 8 


16^ X 32 


16 


0.5846(32) 


0.04531(91) 


0.2659(24) 


5.526(51) 


1.3841(24) 


L = 16 


16^ X 32 


16 


0.537(11) 


0.0554(20) 


0.219(15) 


5.083(67) 


1.431(15) 


L= 16,r > V3 


8^ X 32 


24 


0.636(12) 


0.0280(33) 


0.2997(88) 


6.94(39) 


1.3503(88) 


L = 8 


16^ X 32 


24 


0.579(40) 


0.065(13) 


0.272(27) 


4.60(41) 


1.378(27) 


L = 16 


X 32 


32 


0.609(12) 


0.0369(36) 


0.2809(90) 


6.09(28) 


1.369(90) 


L = 8 


16^ X 32 


32 


0.611(43) 


0.059(13) 


0.295(29) 


4.79(50) 


1.355(29) 


L = 16 


8''' X 32 


18 


0.618(15) 


0.0288(41) 


0.309(11) 


0.83(19) 


1.311(11) 


L = 8 



TABLE VIL Gluonic observables obtained from the static potential for /3 = 2.4. We note that 
ro/a shows significant volume dependence for Lg = 32. 
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p 


Ls > 24 


Ls > 16 


2.3 
2.4 


0.00086(17) 


0.00026(25) 
0.00098(13) 



TABLE VIII: Fit results for the chiral extrapolation of the gaugino condensate, depending upon 
the range of Lg values used. For /? = 2.3, the quality of the Lg > 16 fit is very poor due to nonlinear 
dependence on nij-es that enters at Lg = 24, as can also be seen from Fig. [3j 



IV. EXTRAPOLATION OF THE GAUGINO CONDENSATE 



One important question is the size of Lg necessary to get into the linear regime where 

(AA) ^ Co + cim^es (4.1) 

is a good approximation. Obviously, this serves as an indicator of where we need to be in 
order to have SYM well-approximated. Thus, the measurement of (AA) vs. mres is an impor- 
tant benchmark for determining the regime in which other SYM phenomena can be studied 
with the DWF lattice approach. Another question is the extent to which co,i are sensitive 
to finite spacetime volume (V4 = V x T in our notation). In fact, we find, interestingly, 
that most of the volume dependence is absorbed into mres- AH of this is clearly seen from 
Figs. [3]and|H One sees that to a good approximation, the 8^ x 32 and 16^ x 32 lattice data 
lie on the same line. The smaller value of rrires on the smaller lattice is most likely due to a 
smaller density of near-zero modes. The chiral extrapolation {rrij-cs 0) of (AA)a^ obtained 
from the fit is given in Table IVIIII A feel for the sensitivity to the fitted range of Lg can be 
seen from the two results we provide for (3 = 2.3, which differ by the minimum Lg that was 
included. In fact, the quality of the Lg > 16 fit is very poor due to nonhnear dependence on 
mres that enters at Lg = 24, as can also be seen from Fig. [31 

According to Table IVIH the value of the lattice spacing a is smaller on the /? = 2.4 
lattice. Thus it is surprising that the extrapolated value of (AA)a'^ is larger for j3 = 2.4 
than for /3 = 2.3. On the other hand, we note from Table UTTl that the f3 = 2.4 condensates 
measured in physical units ((AA)rQ) are significantly larger than the (3 = 2.3 ones, given in 
Table [TTl A plausible interpretion is that there are larger renormalizations of the condensate 
on a finer lattice (^ = 2.4), a hypothesis that we are preparing to test with nonperturbative 



renormalization 



42 



43| in an upcoming study. Under this interpretation we have an under- 
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FIG. 3: Condensate vs. rrires for (3 = 2.3, in bare lattice units. Dashed lines show the two linear fits 
(differing by the minimum Lg included) . Extrapolated values together with fit errors are shown at 

m-res = 0. 
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FIG. 4: Condensate vs. mres for = 2.4, in bare lattice units. The dashed line shows the linear 
fit. The extrapolated value together with fit error is shown at mres = 0. 
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standing of why in Table IVlnl the chirally extrapolated value of the lattice units condensate 
for (3 = 2.4 is larger than the one for j3 = 2.3, which does not agree with naive expectations. 
It can also be seen, for instance from the /3 = 2.4, Lg = 32 data in Table IIIIl that a very 
significant finite size effect occurs for the condensate expressed in physical units, a reflection 
of the flnite size dependence of the static potential flt results for this choice of parameters 
(cf. Table rvTTjl . Finally, it is amusing that the extrapolated value of the f3 = 2.3, Ls > 24 
flt is within Icr of the extrapolated value of the (3 = 2.4 flt, though there is no reason why 
this should be true. 



V. CONCLUSIONS 

We have performed a detailed Monte Carlo simulation study of the supersymmetric limit 
of A/" = 1, SU(2) super- Yang-Mills using a DWF lattice formulation. Our work follows from 
an earlier calculation by FKV but signiflcantly extends that work; we use larger lattices with 
two lattice spacings and are able to probe much closer to the chiral limit. Our results for 
the gaugino condensate show the correct theoretical dependence on the residual mass and 
allow for a reliable extrapolation to the chiral limit. Our results provide strong evidence 
for a nonzero gaugino condensate in the supersymmetric continuum limit, and establish 
important benchmarks for future studies. 

Future work that is envisioned is aimed at developing a deeper understanding of the 
conflgurations that are responsible for generating the nonzero gaugino condensate. In par- 
ticular, we would like to elucidate the continuum picture on the cylinder x 5*^, where it 



is monopoles and "KK monopoles" that combine to yield the inflnite volume value jl2|. 

At the same time, two important studies need to be done in order to further develop 
the lattice results presented here. First, we will make a more accurate determination of the 
Sommer scale ro, using the force technique. In that study, smearing of the gauge links and 



other reflnements appearing in 



38| will be used. Second, nonperturbative renormalization 



of the gaugino condensate will be performed [42, l43 |. 



Finally, we note that rather large Lg values were required in order to get mi.es''"o ~ 1/4. To 
improve the situation we envision switching to simulations with modifled versions of DWF 



that have superior chiral behavior 



44 



45|. 
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APPENDIX A: LATTICE FORMULATION 

In this appendix, the A/" = 1, SU(2) SYM lattice action and operators associated with 
;he gaugino condensate are described. The DWF formulation for this theory is identical to 



23|, 



24 



3J], and is described here for completeness. The lattice consists of an SU(2) gauge 
theory with a single Majorana fermion in the adjoint representation. As such, the fermionic 
part of the path integral is the analytic square root of the Dirac determinant, yielding the 
Pfaffian of the corresponding matrix. 
The partition function is 

Z = j [dU] J [d^] J [d^]e-^. (Al) 

Here U^{x), fi = 1,2,3,4 is the four-dimensional gauge link field in the fundamental repre- 
sentation, s) is a (real) five-dimensional Majorana spinor in the adjoint representation 
and $(x, s) is a (real) five- dimensional bosonic Pauli Villars (PV) field with the same indices 
as the Majorana fermion. The coordinate x denotes sites in the four-dimensional spacetime 
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box, which has L sites in spatial directions and T sites in time. The boundary conditions 
along these directions are taken to be periodic for the gauge link fields but antiperiodic 
in time for the fermion and Pauli-Villars fields. The coordinate of the fifth direction is 
s = 0, 1, . . . , — 1, where is the size of that direction and is taken to be an even number. 
The action S is given by: 



S = SciU) + U) + Spvi^, U) . 



(A2) 



SciU) is the pure gauge part and is defined using the standard single plaquette action of 
Wilson: 

5G = /3j](l-^ReTr[f/,]) (A3) 

p 

where j3 = A/g^ and g is the gauge coupling. 
The fermion part Sf{^ ■, U) is given by: 



Sf = - ^{x,s)Df{x,s-x',s')^{x',s') 



where Dp is the DWF Dirac operator in the form of 46| : 



Df{x, s; x\ s') = 5{s - s')p{x, x') + p^{s, s')5{x - x') 



(A4) 



(A5) 



P{x, x') = IY1 + + A - X ) + (1 - i,)Vi{x')S{x' + A - 

+ (mo - 4)(5(a; - x') 



X 



(A6) 



Pndil - s') - mfPiSiLs -1-s')- S{0 - s') s = 

Pr6{s + 1-s') + Pl6{s - 6{s - s') < s < - 1 (A7) 

-mfPR6{0 - s') + PlS{Ls -2-s')- 6{L, -1-s') s = L, - 1 

1±75 



R,L 



(A8) 



where V is the gauge field in the adjoint representation. It is related to the field in the 
fundamental representation: 



K{x)U = lTT[Ul{x)a'^U,{x)a'']. 



(A9) 



with a" the Pauli matrices. In the above equations mo is a five-dimensional mass representing 
the "height" of the domain wall and it controls the number of light fiavors in the theory. In 
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order to get one light species in the free theory one must set < mo < 2 [2l|] . The parameter 
TTif exphcitly mixes the two chirahties and as a result it controls the bare fermion mass of 
the four- dimensional effective theory. In our simulations we have set = 0. 

The fermion field ^ is not independent but is related to \1' by the equivalent of the 



Majorana condition for this 5-dimensional theory 



(AlO) 



where is a reflection operator along the fifth direction and C the charge conjugation 
operator in Eucledean space which can be set to: 



C = 7o72 • 



Therefore, the fermion action can also be written as: 



(All) 



Sf = - '^^{x,s)Mf{x,s;x',s')'${x',s') 



x,x' ,s,s' 



where 



Mp{x, s; x', s') = CR^Dp^x, s; x', s' 



(A12) 



(A13) 



is an antisymmetric matrix as can be easily checked 
gives the anticipated Pfaffian: 



23| . As a result the fermionic integral 



(AM) 



Because det^CR^) = 1 one also has that det(Mir) = det^Dp) and therefore: 

Pf (Mir) = ^/detiDp) . 



(A15) 



The Pauli-Villars (PV) action Spy is designed to cancel the contribution of the heavy 
fermions j47|. Viewing the extra dimension as an internal flavor space 47 1 one can see that 
there are Lg — 1 heavy fermions with masses near the cutoff and one light fermion. The 
PV subtraction removes effects of the Lg hea vy p articles in such a way that the overlap 



determinant is obtained in the L, — > oo limit 



23j. The PV subtraction used here is as in 



48( 1 and is given by: 



Spv = $^(x, s)Mp[mf = l](x, s; x', s')^{x', s' 



(A16) 



x,x' ,s,s' 
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The integral over the PV fields results in: 

r^$ig-Spv = i ^ _ (A17) 

The gaugino condensate is measured using four-dimensional fermion fields that are a 
projection of the five- dimensional DWF 46| : 



Ph^(x,0) + Pl^(x,L,-1) 



(A18) 



In the Ls — oo limit of the theory these operators directly correspond to insertions in the 
overlap of appropriate creation and annihilation operators 47 1. 

Using Eq. (lAlOp and (lAlSP the Majorana condition on the four-dimensional fermion field 

is: 

X = X^C . (A19) 

Because this is the correct condition for a four- dimensional field one can see that the def- 
inition in Eq. (lAlOl) not only produces an antisymmetric fermion matrix Mp but is also 
consistent with the projection prescription in Eq. flAlSp as expected. 



APPENDIX B: RESIDUAL MASS 



Residual chiral symmetry breaking is understood through the axial Ward identity on the 
DWF lattice [s^, which we write here for the case of vanishing bare mass: 

^Ml{^)P'm = 2{JUx)P\Q)) (Bl) 

where is a DWF version of the axial current with isospin index a = 1,2,3 and P"^ is 
a DWF version of the corresponding pseudoscalar current, the interpolating operator for 
the "pion." We note that in the present case the pion is in the adjoint representation of 
the SU(2) gauge group (not to be confused with isospin). The pseudoscalar current J^g is 
different in that it represents a pion at the middle of the fifth dimension: 

J^gix) = -i^ix, Ls/2 - l)PLcrXx, Ls/2) + ^(x, L,/2)PijaXa;, LJ2 - 1). (B2) 

It accounts for the difference in how the left- and right-handed "quark" fields in the DWF 
description transform, which causes a mismatch midway between the domain walls and hence 
the explicit nonconservation of the axial current in the lattice theory at finite Lg. 
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Since in the continuum limit the Ward identity must transition to the continuum form, 
we see that J^q must be related to the pseudoscalar current P through the residual mass 
that is a consequence of 7^ oo: 

J^^ ^ m,,,P\ (B3) 

Thus to extract the residual chiral symmetry breaking one studies the large (imaginary) 
time behavior of the ratio: 

ruj-es = hm = ; — — z — — . B4 

In all of our work we find that this quantity reaches a plateau in the range 3 < t < T — 2, 
and fit to a constant in that region of t. 

APPENDIX C: SIMULATION 



Here we make a few brief remarks on the computational aspects of this 



pro 



ect. Config- 



urations were generated with the rational hybrid Monte Carlo algorithm |49l . 150, iMfl. All 
simulations were performed on Rensselaer's Computational Center for Nanotechnology Inno- 
vations cluster of 16 IBM BlueGene/L machines. We typically used the full capacity of two 
such machines 24 hours/day, and generated configurations for approximately ten months, 
for a total of approximately 30 Million IBM BlueGene/L core hours. The time required 
for data analysis was a small fraction of this, by comparison. Naturally, the large lattices 
(16'^ X 32) with small mj-es values were the most costly to generate. 

In the rational approximation used to generate configurations, we found that it was 
necessary to go to rather high degrees, due to a very wide spread between lowest and highest 
eigenvalues of the Dirac operator. This occured because we performed our simulations at 
vanishing bare fermion mass mj = 0, relying on the finite but large value of Lg as an 
infrared regulator. Typically, the Metropolis step required degrees between 15 and 20 in 
the computation of the change in the Hamiltonian, in order to have sufficient accuracy. 
Moreover, it was not unusual to require between 50 and 100 steps in the leapfrog integration 
for a trajectory of r = 0.5 simulation time units, in order to get reasonable acceptance 
rates at large Lg. Again, this was a result of small eigenvalues of the Dirac operator. 
Naturally, these features led to very slow updating. For the Lg values that we simulated, a 
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/? 


y X T 


Ls 




configs./day 


2.3 


8^ X 32 


48 


60 


140 


2.3 


8^ X 32 


64 


80 


70 


2.3 


16^ X 32 


32 


90 


40 


2.4 


8^ X 32 


48 


100 


80 


2.4 


16^ X 32 


48 


200 


11 



TABLE IX: Example timing results on a single BlueGene/L rack (1024 dual core nodes). Here, 
Nr is the number of steps in the leapfrog trajectory. 

single BIueGene/L rack was able to produce (9(10) to (9(100) configurations per day. Some 
examples are given in Table IIX[ The last row represents a run that was not reported in 
the main text as it was too slow for a reasonable data set to be generated in a practical 
time-frame. As a result, our forthcoming work will set rrif ^ and perform an ^ 
extrapolation for these larger volume, large simulations. 
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